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Abstract 

We show that an expansion of nuclear forces about the chiral hmit is formally 
consistent and is equivalent to KSW power counting in the ^Sq channel and Weinberg 
power counting in the ^Si — ^Di coupled channels. Numerical evidence suggests that 
this expansion converges. The feasibility of making contact between nuclear physics 
and lattice-QCD simulations is discussed. 



1 Introduction 



The last decade has seen significant progress toward the formulation of an effective field 
theory (EFT) description of multi-nucleon systems (for recent reviews see Ref. |jl|, ^, 
^ ^ ^). The importance of uncovering such an EFT cannot be overstated as it is 
this theory alone that will allow for rigorous calculations of both elastic and inelastic 
processes in multi-nucleon systems in a framework consistent with the Standard Model 
of strong and electroweak interactions. An inherent advantage of an EFT framework 
is that the uncertainty associated with the computation of any given observable can be 
estimated and controlled with the power-counting scheme that defines the theory (along 
with the regulator). The fundamental difficulty in formulating power-counting schemes 
for multi-nucleon systems is that the interaction of two or more nucleons near threshold 
is intrinsically nonperturbative. This manifests itself through the presence of scattering 
lengths that are much larger than one would naively expect from QCD, ~ AgJ^^ For 
very low-momentum processes, such as np — > dj in the nucleosynthesis region, a consistent 
and converging EFT, EFT(7f), has been developed to describe strong and electroweak 
processes in multi-nucleon systems [0. However, attempts to develop a fully consistent 
theory to describe processes involving momenta larger than the mass of the pion have so 
far been unsuccessful. 

Weinberg's (W) original proposal @ for an EFT describing multi-nucleon systems was 
to determine the nucleon-nucleon (NN) potentials using the organizational principles of 
the well-established EFTs describing the meson-sector and single-nucleon sector (chiral 
perturbation theory), and then to insert these potentials into the Schrodinger equation to 
solve for NN wavef unctions. Observables are computed as matrix elements of operators 
(that have their own chiral expansion in this EFT) between these wavefunctions. W power 
counting has been extensively and successfully developed during the past decade to study 
processes in the few-nucleon systems (see, for example, Ref. 0). This method requires 
numerical solution of the Schrodinger equation and is similar in spirit to traditional nuclear 
physics potential theory. Unfortunately, there are formal inconsistencies in W power 
counting [Q, in particular, divergences that arise at leading order (LO) in the chiral 
expansion cannot be absorbed by the LO operators. Problems persist at all orders in the 
chiral expansion, and the correspondence between divergences and counterterms appears 
to be lost. This formal issue was partially solved by Kaplan, Savage and Wise (KSW) 
who introduced a power counting in which pions are treated perturbatively [Q. The NN 
phase-shifts and mixing angle in the ^5*0 and ^5*1 — ^Di channels have been computed 
at next-to-next-to-leading order (N^LO) in the KSW expansion by Fleming, Mehen and 
Stewart (FMS) ||T0| from which it can be concluded that the KSW expansion converges 
slowly in the ^5*0 channel and does not converge at all in the ^5*1 — ^Di coupled channels. 
FMS found that the part of the scattering amplitude that survives in the chiral limit is 
large enough to destroy convergence in the ^5*1 — ^Di channels. This is in contrast with 
converging numerical results found using W power counting in this channel [T^ . 

Both W and KSW power countings have desirable features that one would hope exist 
in the, yet to be constructed, EFT that describes multi-nucleon systems. It is therefore 
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important to consider two fundamental questions in this area: 



1. Where is W power counting consistent? It has been argued that there are funda- 
mental problems with the chiral expansion and the momentum expansion resulting 
from W power counting |p. In this paper we confirm that W power counting is 
formally inconsistent in the ^5*0 channel, but we do find that W power counting is 
consistent in the — ^Di coupled channels. The arguments presented in Ref. ^ 
were based on a perturbative analysis, whereas the correct renormalization of sin- 
gular potentials is intrinsically nonperturbative |jT3|, 0, ^ . As we will see the same 
is true in the ^Si — ^Di coupled channels. 

2. Where is pion exchange perturbative? The analysis of FMS demonstrated that 
pions, or more explicitly, the Yukawa part of the NN potential and radiation pions, 
when treated perturbatively give rise to a converging expansion for the ^5*0 scattering 
amplitude. They also showed that one-pion-exchange (OPE) in the ^Si—^Di coupled 
channels is not perturbatively convergent. In this paper we confirm the result of 
Ref. [§] that part of the pion contribution to scattering in the ^Sq channel requires 
perturbative treatment in order to obtain a consistent theory. However, there is no 
such requirement in the ^5*1 — ^Di coupled channels. 

In this paper we will argue, supported by numerical evidence, that an expansion about 
the chiral limit provides a systematic power counting for multi-nucleon systems. This 
expansion is found to be equivalent to KSW power counting in the ^5*0 channel and 
equivalent to W power counting in the ^5*1 — ^Di coupled channels, i.e. it selects only 
the desirable features of both power countings. To understand the physics behind this 
expansion, it is sufficient to recall that in the ^Si — ^Di coupled channels, OPE includes 
a strong tensor component at short distances that persists in the chiral limit and is 
large enough so that it must be summed to all orders, as required by W power counting. 
However, deviations from the chiral limit V{r\ m-j^) — V{r; 0) can be treated perturbatively 
in all channels, and in fact, such a perturbative expansion is required in the ^5*0 channel 
but not in the '^5'i — channel. 

In section 2 we discuss the ^5*0 channel. In particular, we review the problems with 
W power counting in this channel. Using the Standard Toy Model of Nuclear Physics 
(the three- Yukawa model developed in Ref. ||16|), detailed in the Appendix, we show that 
the Yukawa part of OPE can be treated perturbatively if the short- distance physics is 
properly accounted for in the EFT. We outline a method for improving the convergence 
of the ^5*0 phase-shift in the realistic case. In section 3 we investigate the ^5*1 — ^Di 
coupled channels. The LO amplitude in W power counting is analytically renormalized 
and a comparison with numerical data is presented. The NN phase-shifts are computed 
at high orders in perturbation theory to study the convergence of the V{r; m^^) — V{r] 0) 
expansion. In section 4 we discuss possible scenarios for the light-quark mass dependence 
of the deuteron binding energy and NN scattering amplitudes, and dimly illuminate the 
long and winding road between nuclear physics and lattice QCD. We conclude in section 
5. 
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2 Scattering in the ^/So-Channel 

NN scattering in the ^5*0 channel has been studied extensively in EFT investigations as 
it results from a relatively simple interaction between nucleons. OPE provides the long- 
distance Yukawa interaction and there is a large short-distance interaction, affectionately 
known as the hard-core. This structure makes the channel a useful testing ground for 
new ideas. 

The ^5*0 channel has been investigated using both W and KSW power countings with 
interesting results. First, calculations up to N^LO in W power counting have been per- 
formed and reproduce the phase-shift data in this channel and in many other 
channels very well for a range of ultra-violet (UV) momentum cut offs. However, there is 
a sensitivity to the momentum cut off due to an inconsistency in W power counting M. 



Second, calculations up to N LO in KSW power counting have been performed & 10], and 



are found to converge to the results of the Nijmegen phase-shift analysis OT^. The con 



vergence of the expansion can be improved substantially by promoting the effective range 
from NLO to LO in KSW power counting |jl6|, implemented via the dibaryon field 1^ 



or the z-parameterization pO |. 



2.1 Problems with W Power Counting in the ^5*0 Channel 

In W power counting the momentum-independent four-nucleon interaction and OPE form 
the LO NN interaction. The Schrodinger equation is used to determine the LO phase- 
shifts for NN scattering from this LO interaction. This procedure corresponds to summing 
an infinite number of insertions of the LO interaction. Naively, this would appear to be 
a sensible procedure, however, subtle problems arise in solving the Schrodinger equation 
due to its bad UV-behavior and it must be renormalized, as discussed in detail in Ref. . 
The OPE central and tensor potentials are 

Vc{r; m^) = -a^ ml (1) 

r 

and 



Vrir; m,) = -a, ml — — ( 1 + — - + ] , (2) 

respectively, where 



m^r rnir"^ 



and we use gA = 1.26, F^^ = 93 MeV and M = 940 MeV. Scattering in the ^5*0 channel 
depends only on the central potential, Vc- Near the origin, Vc is Coulombic [] and gives 
rise to NN wavefunctions that are a linear combination of regular and irregular Bessel 
functions. At LO in W power counting, OPE is renormalized by a single delta-function 

^The 5^^\r) contribution to Vc is not shown explicitly but is absorbed into the definition of the local 
four-nuclcon operator defined in eq. m 
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counterterm, Cq. For the purposes of this paper, we choose to smear the delta-function 
counterterm(s) over a finite volume near the origin and for computational ease choose a 
square well of radius i?, 

C««'"M-^^^||iF^-l'o9(fl-r). (4) 

The OPE potential is similarly cutoff at the radius R = 1/A. The advantage of this 
simple regulator is that analytic matching can be performed 0, and at zero center-of-mass 
(CoM) energy one finds 

^-MVocot (^^-MVoR^ = -mlMa^ log (^-^) + 0{R) (5) 

where i?* is an intrinsic length scale to be determined numerically. Treating v^— MVq-R as 
an expansion in mT^, we find the depth of the energy-independent square well as a function 
of R to be 

V.iR: n) ^ -i2n ^ If^ - H!f - log (|) + O(R'). (6) 
where n labels the branch of the cotangent. The presence of a multi-branch structure is 




R (fm) 

Figure 1: The solid line represents the running of Vq as a function of the cutoff R (in 
fermis), taken from eq. with n = 1, for the physical value of rriT^, while the dotted line 
neglects the R^^ part of the running. The dots are extracted directly from a numerical 
solution of the Schrddinger equation for the physical value of m.^. 

related to the accumulation of bound states inside the square well. Of course the presence 
of unphysical bound states is innocuous as long as the binding energies of such states 

^ The physics we obtain will be regulator independent. Nevertheless, we have been unable to dimen- 
sionally regulate the '^S'l — '^Di coupled channels. 
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are near the cutoff of the EFT. In fig. |I| the R dependence of Vq, as given by eq. (Q), is 
compared to points extracted from the numerical solution of the Schrodinger equation. 
For a given i?, Vq is tuned to give the scattering length, = —23.714 fm, and in fig. |^ 
we show the ^5*0 phase-shifts for different values of R. In momentum-space notation the 
matching equation takes the form: 



Co{^\n) + mlD2{^) = {2n + if 



3MA 



3A2 



log 



A* 

a" 



(7) 



and therefore the D2 operator, which is formally subleading in W power counting, must 
be promoted to LO, in agreement with the perturbative argument of Ref. Although 
the logarithmic divergence is suppressed by a power of A compared to the leading term 
in eq. |^, it is a true divergence in physical quantities which must be renormalized at 
leading order in W power counting (see eq. |TD| below). This offers a proof that in the 
^5*0 channel, an expansion of the action to a given order in W power counting does not 
consistently remove all cutoff dependence to that order P2|. Of course the D2 contribution 
is numerically small, as demonstrated by the dotted curve in fig. which neglects the A"^ 
corrections to the running. This smallness explains why W power counting has been found 
to work well in this channel over a moderate range of cut offs [jll], 0, |^. 




Figure 2: The ^Sq phase-shift plotted versus CoM momentum. The solid lines are the 
Nijmegen phase-shift analysis /|7^. The long-dash lines corresponds to R = 0.2 fm (A = 
985 MeV), the medium-dash lines correspond to R = 0.4 fm (A = 492 MeV), and the 
dotted lines correspond to R = 1.4 fm (A = 140 MeV). The left panel corresponds to 
setting V2 = 0, while the right panel includes V2 in such a way to reproduce the measured 
effective range, r '^'^ . 



How does this problem manifest itself at the level of Feynman diagrams f^? In dimen- 
sional regularization the divergent part of diagram (a) in fig. ^ is 



2567r2F2 
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where C is the coefficient of the four-nucleon operator, and n = 4 — 2e is the number of 
space-time dimensions. This diagram requires a counterterm, D2-, with a single insertion 
of the hght quark mass matrix, m^, to render the amphtude finite, however, such coun- 
terterms occur only at higher orders in the chiral expansion. Thus, the inconsistency of 
W power counting is unmasked diagrammatically. 

(a) (b) 

Figure 3: Graphs with logarithmic divergences at (a) LO in W power counting and NLO 
in KSW counting and (b) NLO in W power counting and N^LO in KSW counting. The 
solid lines are nucleons and the dashed lines are pions. 

It is instructive to understand how the divergence in eq. (P) arises in position space 
with a square-well regulator. The radial wavefunctions resulting from resumming the LO 
bubbles are phase-shifted plane- waves of the form 

u{r) ~ sin [kr + (5lo] , (9) 

where 5lo is the LO phase-shift. Perturbatively inserting the difference between the OPE 
potential for rriq ^ and its value in the chiral limit, and regulating the UV divergence 
with a square well of a radius R much smaller than the scattering length R <^ a, the NLO 
contribution to the scattering amplitude is 

{ij'\V{m^)-Vmi') ~ glml r dr sin^kr + Si^o] ^^^^^^^^ 

J R T 

~ -g\mlsw? 5i^o log (m^i?) + ... , (10) 

where the ellipses denote terms that are finite as -R — 0. At the order at which y(m^) — 
V{0) enters, there must be a contribution from a local operator of the form log (R/R^,), 
where R^, is some intrinsic length scale, in order to render the sum finite and independent 
of the renormalization scale R. 

By the analysis above, one is naturally led to the conclusion that the contribution from 
OPE and the O {rriq) counterterms must contribute to any amplitude at the same order in 
the expansion, and this is what leads to KSW power counting. However, a more general 
conclusion to draw is that the difference between the OPE contribution for rrig ^ and 
the OPE contribution in the chiral limit must occur at the same order as the O iva^ 
counterterms. In most cases, these two conclusions yield identical amplitudes, however, 
in the ^5*1 — ^Dx channel they do not, as we will see subsequently. 

^ If the LO radial wavefunctions behaved differently as r ^ 0, the divergence structure would be 
modified. Such a modification can only arise from a re-ordering of the perturbative expansion, and 
consequently the hierarchy of higher-dimension operators as determined from naive arguments will also 
be modified. 
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2.2 Are Pions Perturbative? 



In KSW power counting the renormalization group (RG) flow of tlie LO four-nucleon in- 
teraction promotes it to one lower order in the chiral expansion, while the OPE interaction 
does not evolve. The LO NN scattering amplitude arises from resumming the bubbles 
with an infinite number of insertions of the four-nucleon interaction. At NLO, there are 
contributions from OPE, insertions of O {rUq) four-nucleon operators and from insertions 
of (9 (|pp) operators. Operationally, this hierarchy arises from the assignment that the 
scattering length scales as an inverse power of the small expansion parameter, a ~ l/Q, 
while all other parameters are natural size. The scattering amplitude in the ^5*0 channel 
has been computed to N^LO in KSW power counting in Ref. [1^, where it was found to 



converge slowly. However, the cause of the slow convergence with KSW in this channel 
stems from the treatment of the short- distance contribution to the scattering amplitude 
and not from perturbative OPE. All momentum dependence arising from the hard-core 
of the NN interaction is treated in perturbation theory and as these effects are large, even 
over a relatively small momentum range, the convergence is slow. 

2.3 Lessons from a Model 

A toy model that has proved useful in understanding the structure of the NN interaction is 
the three- Yukawa model 0. One includes a Yukawa pion, rho and sigma meson, with their 
physical masses (the sigma meson is given a mass of 550 MeV, despite its non-existence in 
nature) and their couplings are adjusted to reproduce the scattering length and effective 
range in the ^5*0 channel. Convergence problems with KSW power counting were high- 
lighted with this model by Steele and Kaplan [|T6l , and it is important to understand the 



cause of these problems before proceeding further. 

Several valuable lessons about both KSW and W power counting formulations can be 
learned from the STMoNP. First, pions are perturbative in the ^So channel. The full 
theory defined via the Schrodinger equation can be systematically perturbed in the num- 
ber of insertions of the pion potential. In addition, one can recover these perturbative 
amplitudes with an EFT construction provided one retains a sufficient number of countert- 
erms. Second, the fact that W power counting | converges nicely, while KSW converges 
quite slowly [|^, and appears to encounter problems at a scale ~ m-,, is now understood. 



The problem is not with the perturbative pions component of KSW but rather with the 
treatment of the short- distance physics. By resumming the "hard-core" into a dibaryon 
field, this modified-KSW (or modified- W) power counting does provide a perturbatively- 
convergent amplitude. The peculiar feature, however, is that by only including a small 
number of local counterterms and by eyeballing the phase-shift plots, one would guess 
that the breakdown scale for the theory is ~ m^r. However, a more systematic analysis. 



"^We will refer to this model, which we detail in the Appendix, as the Standard Toy Model of Nuclear 
Physics (STMoNP). 

^In the STMoNP, W power counting does not lead to any ambiguities because the long-distance 
component of the nucleon interaction is sufficiently UV convergent. 
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such as that provided by error plots or by going to higher orders in the expansion, reveals 
a breakdown scale of mo-/2, set by the radius of convergence of the effective range ex- 
pansion for the short-distance physics. Third, it is clear from this analysis that W power 
counting is overkill. One does not need to resum pions exchanges in order to arrive at a 
meaningful phase-shift for this theory, and in fact, closed-form analytic expressions exist 
and reproduce the numerical solutions of the Schrodinger equation up to rria/l. 

Note that the agreement with data (in the left panel of fig. breaks down at very-low 
momenta because the scattering length is large and the effective range in this channel is 
large as well, r "^o = 2.73 fm, comparable to twice the pion Compton wavelength. This 
suggests that the two-derivative C2 operator should be promoted to LO, and is easily 
included in our analysis via an energy-dependent contribution to the depth of the square 
well: 

C2VH^^\r) C2ME8^^\r) ^ ~ = fcVa e{R - r) (11) 

where E = k^/M is the CoM energy. The first step amounts to a field redefinition which 



removes momentum dependence in favor of energy dependence [T^, |2^, M. The value of 
V2 as a function of R, is chosen to recover the experimental value of r^". In the right 
panel of fig. ^ we see the much-improved phase-shift with resummed range corrections, 
providing an example of the manner in which the momentum expansion can be tailored 



to experimental realities [jT9|. 



3 Scattering in the ^Si — ^Di Coupled Channels 

In addition to the long-distance Yukawa OPE interaction and the hard-core interaction, 
NN scattering in the '^5'i — ^Di coupled channels has a strong tensor component. At dis- 
tances r -C m-jr the central potential is negligible and the nucleons interact entirely through 
the singular tensor interaction. The chiral expansion of the tensor force is equivalent to 
the short- distance expansion and is given by 



VT{r;m^) = _^ + !^ + o{mt) . (12) 
r'^ 2r 



It is the leading term in this expansion which leads to the breakdown of the KSW expan- 
sion in the ^5*1 — ^Di coupled channels [|T^]. Note that a significant consequence of chiral 
symmetry is that there is no 1/r^ term in the short-distance expansion of Vrir; m.,^). One 
may hope that the convergence problems in the ^Si — ^Di coupled channels can be cured 
as in the ^5*0 channel by a re-ordering of local operators representing this short distance 
physics. Unfortunately, this is not the case. For potentials of the form l/r", with n < 3, 



k cot 6 is non-analytic in k (for a comprehensive review of singular potentials see Ref. |p,4|| ) 
and consequently the chiral limit of Vrir; 771^^) cannot be represented by local operators, 
and must be summed to all orders. At first glance, this might seem to pose an insur- 
mountable problem for the EFT effort. KSW have shown that if one isolates terms in 
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perturbation theory with insertions of the tensor force, it appears that an infinite number 
of counterterms are necessary to properly renormahze the highly- singular behavior. This 
was motivation for treating the pion perturbatively. Fortunately, this argument is not 
correct, essentially because it relies on perturbation theory, and the correct renormaliza- 
tion is invisible in perturbation theory. This is similar to what occurs in the three-body 
system with short-ranged interactions |^ . In recent work, some of us showed that a 1/r" 
attractive singular potential can be renormalized by a single energy- independent square 



well |T^. Here we will extend this result to the realistic case, which amounts to show- 
ing that the momentum expansion assumed by W power counting is consistent in the 
^5*1 — ^Di coupled channels. The investigation is in similar spirit to the work in Ref. ||26[| 
who consider an OPE potential and work with a subtracted T-matrix, and to the work 
in Ref. which examines scattering and electromagnetic observables in the ^5*1 — ^Di 



coupled channels. 

3.1 W Power Counting at LO 

At LO in W power counting, OPE in the ^Si — ^Di coupled channels is renormalized 
by a single delta-function counterterm, as it is in the ^S'o-channel. For demonstrative 
purposes, we again choose to regulate the delta-function with a square well, which has 
been implemented previously in very forward-thinking work by Sprung et al The 



S'-wave and D-wave components of the J = 1 system, u{r) and w{r) respectively, are 
coupled by the tensor potential, Vr, defined in eq. (0). The potential outside the square 
well is 

' 1^ -2^/2 MVrir; m^) -M (Vb(r; m,) - 2VT{r; m,)) - G/r^ J ' ^ ""^ 

while the potential inside the square well is 

, . / -M{Vo + eV2) \ 

^^^"^=1, -MiVo + k'V,)-6/r^ ) ■ ^^4) 

Here we have included V2, but of course it is not necessary at LO in W counting. We have 
also left out ^Si — ^Di superscripts which differentiate this square well from that which 
appears in the ^Sq channel. Defining \1/ to be 



the regulated Schrodinger equation is 

^"(r) + (A;2 + Viir; m^)e{r - R) + Vs{r)e{R - r)) ^ = 0. (16) 

Renormalization concerns the behavior of physics at short-distances, or more specifically, 
is due to our lack of understanding of physics at short- distances. If i? < r <C Ma^ = 
1/Ann, then we can neglect the angular- momentum barrier and keep only the chiral 
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limit of the tensor force. Moreover, for \k\ <^ ^nn the total energy can be treated 
as a perturbation. In this short-distance limit, VL(r;m^) can be diagonalized and the 
Schrodinger equation can be solved exactly. In the diagonal basis the Schrodinger equation 
decouples into an attractive singular potential and a repulsive potential. The solution for 



the attractive singular potential, u, is a linear combination of Bessel functions |T5 



and it is useful to think of the short- distance limit of the Bessel functions using the WKB 



approximation [|T^. The first two orders in the WKB approximation reproduce the LO 
asymptotic expansion of the Bessel functions, while NLO in WKB gives the leading energy 
corrections. The wavefunction at this order is 



k) ^ A .'/'cos {2f-^ - + + .v.] . (17) 

where A is a dimensionful normalization constant and 0o and 0i are the coefficients of 
a Taylor-series expansion of the asymptotic phase Note that this wavefunction 

vanishes at the origin. 

The solution for r < i? is simply that of a square well of height Vq + k'^V2. Matching 
logarithmic derivatives of the interior square well and exterior attractive solutions at 
r = R yields the matching equation 




-MVoRcot { \/-MVoR ) = 7 + ^/^^tan hJ^^ + 



for Vn and 



2 



MVo 



~MVo cot J-MVqR + MVqRcsc' J-MVqR 




tan 2W— -^ + 00 - TR+-f^\—fT^h sec^ 2W— -^ + 0o (19) 



for V2. These equations determine the renormalization-group flow of Vq and V2. Note that 
this flow is multibranched and non-analytic in the coupling f\ The phases 0o and (pi 
are determined numerically. For a given R, Vq is tuned to recover the physical value of 
the scattering length in the ^5*1 channel, "^^^ = 5.425 fm, from which eq. ( [T8| ) determines 
the phase 0o- ^2 is chosen to recover the correct physical value of the effective range, 
^(^■Si) _ X 749 fj-Q from which eq. ( ]T9| ) determines 0i. In flg. ^ the square-well depths Vq 
and V2 are shown as functions of the square- well width R taken from eq. (p!8[) and eq. (plOl). 
The LO phase-shifts in the ^5*1 — ^Di coupled channels are shown in flg. |^ and flg. ^ with 
V2 = 0. While the phase-shifts 6q and 62 are clearly independent of the square- well radius 
R, the mixing parameter ei exhibits R dependence. However, an error plot of ei indicates 
that the R dependence and the deviations from the Nijmegen phase-shift analysis are 
higher order in the momentum expansion, beginning at O (/c^) (and do not indicate an 
inconsistency in the power counting). 
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Figure 4: The solid lines represent the running of MVq and MV2 as a function of the 
cutoff R, taken from eq. ( |7^ and eq. i^TV{), respectively. The dots are extracted directly 
from a numerical solution of the Schrddinger equation. The different branches in the left 
panel correspond to a different number of nodes in the square well, i. e. an artifact of this 
particular regulator. The curves continue to smaller values of R, but we have not shown 
them. Further, analogous branches exist for MV2, but we have not shown them. For each 
branch, the ultraviolet phase is fit to the smallest R data point to produce the theoretical 
curve. 



I P I (MeV) 




Figure 5: The ^Si — ^Di phase-shifts as a function of the CoM momentum. The solid line 
is the Nijmegen phase-shift analysis The long-dash line corresponds to R = 0.45 fm 
(A = 438 MeV), the medium-dash line corresponds to R = 0.21 fm (A = 938 MeV), and 
the dotted line corresponds to R = 0.10 fm (A = 1970 MeV^. The phase-shifts 60 and 62 
are clearly cutoff-independent (see text). 



The scaling of the leading explicit chiral symmetry breaking operators in the UV can 
be found by perturbatively inserting the difference between the OPE potential for rriq 7^ 

^For sufficiently small coupling, analytic behavior in 0;,^ is recovered, as detailed in Ref. [p5[. 
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Figure 6: The left panel shows ei as a function of the CoM momentum.. The solid line 
is the Nijmegen phase-shift analysis [P^J . The long-dash line corresponds to R = 0.45 fm 
(A = 438 MeV), the medium-dash line corresponds to R = 0.21 fm (A = 938 MeV), and 
the dotted line corresponds to R = 0.10 fm (A = 1970 MeV). The right panel shows an 
error plot for ei. Momentum \p\ is in units of MeV . The quantity A is the difference 
between the EFT calculation discussed in the text and the Nijmegen phase-shift analysis. 
The slopes of the curves for R = 0.25 fm and R = 0.45 fm are ~ 2. 



and its value in the chiral limit, 



Ul)'\V{m^) - Vmi)) ~ glA^ / dr cos^ 

Jr 

where the ellipses denote terms that are finite as i? — >■ 0. From the momentum-space 
perspective these leading explicit-breaking effects have a strange A~^/^ scaling where A^^, 
is the chiral symmetry breaking scale. In practice, R is kept finite and while one might 
consider putting in a D2 operator in order to remove this R dependence, it is formally not 
necessary. This is in contrast to the ^5*0 channel, where logarithmic dependence is found, 
eq. (p!0[) , which does require the inclusion of a D2 operator. 

3.2 - ^Di Phase-Shifts Perturbatively in V{r] m^) - V{r; 0) 

There is now reason to believe that a perturbative expansion about the chiral limit is 
formally consistent in both the ^Sq and ^5*1 — ^Di channels. In order for this expansion to 
be useful it must also converge. At NLO there will be finite 0{m\) corrections from the 
perturbative expansion in V^(r; m^r) — V(r;0). However, since is itself an expansion 
about the chiral limit, the D2 operator will contribute distinctly from the Co operator. 
Hence, knowledge of D2 requires a distinct measurement involving multi-pion interactions 
with nucleons. In the absence of knowledge about D2, a well-defined, and significant com- 
putation that can be performed is the perturbative expansion of the scattering amplitudes 



~ -g'^A'miR-^ + 



,(20) 
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in V{r;mT^) — V{r;0), the difference between the OPE potential evaluated at the phys- 
ical value of and the OPE potential for a massless pion. Unfortunately, this does 
not correspond to the chiral limit, as the couplings such as qa and are held fixed at 
their physical values [J The STMoNP suggests that one should expect the perturbative 
phase-shifts to be most deformed at lower-momenta, and approach the observed value at 
higher momenta. Fig. |^ shows the phase-shifts and mixing angle of the ^Si — ^Di coupled 
channels order- by-order in perturbation theory out to N^LO. The expansion is conver- 
gent I, but only slowly. Despite the relatively poor convergence, this result is important 




Figure 7: The 60, 62 phase-shifts and ei mixing angle in perturbation theory out to N^LO 
in V{r; rrij,) — V{r; 0) . The solid line is the result of the Nijmegen phase-shift analysis / (7^/ , 
while the dotted curves correspond to the LO, ISPLO, N^LO and N^LO results, as indicated 
on each panel. We have not shown odd orders to reduce clutter. These results were 
obtained for a square well of radius of R = 0.25 fm, with an energy-independent square- 
well depth of MVq = 1.01 x 10^ MeV^ and MV2P = 23.5AP. 

as it demonstrates that the expansion of the potential about the chiral limit is convergent. 
Clearly, one must be concerned about the slow convergence of this expansion. One hopes 
that this will not seriously impact the construction of a systematic chiral expansion, but 
unfortunately it may. 

4 How Nuclear Forces Depend on niq 

One of the important goals of the EFT investigations is to obtain the dependence of nuclear 
observables upon the masses of the up and down quarks, and eventually the strange quark. 
This would allow one to determine the deuteron binding energy in the chiral limit, and 
other interesting but perhaps somewhat academic features of chiral-limit nuclear physics. 
From a practical point of view, knowledge of the quark-mass dependence of observables 
will allow one to isolate terms in the chiral Lagrangian that contribute differently to elastic 

^ An expansion in V{r-, m-,^) (not subtracting the V{r\Q)) should recover the FMS analysis of KSW 
power counting at N^LO ||10| minus the contributions from radiation pions. 

*By convergent, we mean that the perturbative expansion converges to the solution to the Schrodinger 
equation for the full potential V{r] m-^). 



13 



and inelastic processes involving pions. This is an important step toward developing a 
unified theory of elastic and inelastic processes. Further, knowledge of the quark-mass 
dependence of these observables will provide the bridge between proposed lattice-gauge 
theory calculations, performed at quark masses much larger than the physical values, and 
nature f\ 

In multi-nucleon sector calculations, quark-mass dependences are both explicit and 
implicit. For instance, in the NN-sector there is explicit m^-dependence in the mass of 
the pion in the nuclear potential and also in four-nucleon operators. However, there is 
implicit m^-dependence in the couplings in the meson sector (e.g. F^) and single-nucleon 
sector (e.g. qa)- Without measuring processes involving either initial or final-state pions, 
the coefficient of each four-nucleon operator that has explicit rrig dependence cannot be 
determined. Therefore, what follows is for illustration only, and we do not claim to 
presently know the deuteron binding energy, or any other finely-tuned observable in the 
chiral limit. 



4.1 Explicit m^-Dependence only 

It is interesting to know the behavior of observables in nuclear systems as a function of 
the explicit factors of rrig that appear in the nuclear potential. That is to say, the value 
of the pion mass in the OPE potential scales as ~ but other quantities, such as 

Ftt, M and qa are fixed to their measured values. It is this scenario that can be directly 
compared with all potential-model calculations, where the implicit m^-dependence of the 
short- distance part of the nuclear interaction is unknown. 

As we demonstrated earlier, with square-well regularization of the ^5*1 — ^Di coupled 
channels, the m^-dependent four-nucleon counterterm required to yield a renormalization- 
scale (R) independent result scales like D2 ~ R^^'^ as i? — 0, and therefore is naively 
higher order in the power counting. Hence, we may vary the value of the explicit factors 
of in the nuclear potential without an accompanying modification to the depth of the 
regulating square well. A plot of the deuteron binding energy as a function of resulting 
from numerically solving the Schrodinger equation is shown in fig. For the physical 
value of the pion mass we recover the deuteron binding energy to reasonable accuracy, 
= 2.211 MeV (essentially independent of R). In the chiral limit one finds that the 
deuteron is bound by B'^ ~ 4.2 MeV. This value is still somewhat small compared to 
F^/{2M) ~ 10 MeV, that one might expect to arise in QCD, and therefore one would 
conclude that the deuteron is still weakly bound in the chiral limit! Our calculation 
of the deuteron binding energy in the chiral limit agrees with that obtained with the 
AV18 potential of B^/"^^^^' ~ 4.1 MeV. These very-similar results are in qualitative 



agreement (but of opposite sign) with a previous estimate [^3] based on a first-order 
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Recently, remarkable progress has been made in understanding how to compute higher-order coef- 



ficients of the chiral Lagrangian describing the pseudo-Goldstone bosons 1 29 . Partially-quenched lattice 
calculations have been matched onto a partially-quenched chiral Lagrangian and combinations of the 
Gasser-Leutwyler coefficients [Q, the L^'s, have been determined. At this point in time, a partially- 
quenched EFT describing the single and multi-nucleon sectors does not exist. 
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Figure 8: The deuteron binding energy as a function of the m^r that explicitly appears in 
the OPE potential. For R = 0.25 fm, an energy-independent square well depth of MVq = 
1.01 X 10^ MeV^ and a square well that depends linearly on energy of MV2k'^ = 23.54fc^ 
reproduce the scattering length and effective range in the ^Si channel for m^^ = 139 MeV. 



perturbative calculation with wavefunctions from the Paris potential and Bonn potential. 

In addition to the deuteron binding energy, one gains useful insight into the chiral limit 
of the NN-system by examining the phase-shifts and mixing parameters of the ^5*1 — ^Di 
coupled channels. We have considered the perturbative expansion of the phase-shifts and 
mixing angle previously, but it is worth restating in this discussion. By numerically solving 




Figure 9: The phase-shifts 5o, 82 and the mixing parameter e\ of the "^Sx — ^Di coupled 
channels, generated by the parameter set as fig. The solid curves correspond to the 
Nijmegen phase-shift analysis j{P^J, while the two sets of dotted curves correspond to our 
calculated result for the physical value of the pion mass, denoted by m-^ in the plots, and 
the chiral limit, as denoted in the plots. 

the Schrodinger equation we find the phase-shifts and mixing parameter shown in fig. ^, 
see also fig. |^. It is clear that the observed ^5*1 phase-shift does not differ significantly 
from that resulting from taking ttItt to zero in the potential. However, as expected, the 
^Di phase-shift and mixing parameter ei differ significantly from this "chiral limit" . As 
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the deuteron is bound largely by the strong attractive tensor force, the observed increase 
in El explains the deuteron binding energy computed above. 



4.2 Implicit and Explicit mg-Dependence 

In addition to the explicit dependence of the OPE, there is implicit mg- dependence, as 
the couplings in the meson and single-nucleon sector also depend upon at subleading 
orders. Such dependences have been computed previously |^ for QCD with SU{2)i ® 
SU{2)r chiral symmetry. The physical values of the pion decay constant Fj^, the nucleon 
mass M and the axial coupling constant qa are given in terms of the pion mass at lowest 
order in the chiral expansion m, and F, Mq and g, their respective values in the chiral 
limit 



F^ 
M 

9a 



1 + 



m 



U + 0(m 



167r2F2 

Mo - Am'^ci + O (m^) 
2^2 + 1 



^^rn^log(^J+OK) 



(21) 



where h = 4.4 ± 0.2 ||, H, ci 1 GeV~^ H, and we have use m = 140 MeV. We 



have retained only the leading chiral-logarithmic contribution to qa, and have chosen a 
renormalization scale of A = 500 MeV. In addition to these relations, there is also rriq- 
dependence in the relation between the tcNN coupling and the axial matrix element 
computed at LO in Ref. |p6 



as 



g-n-NN 



Ftt \ 9a , 



where we will use rfig ~ —1.0 GeV~ ||3g]. Even with this set of mg-dependences included 
into the calculation, there very well could remain uncalculated rrtg-dependences. In partic- 
ular, the implicit rrig-dependences modify the strength of the chiral-limit OPE potential. 
This will give rise to a stronger R-dependence in the local four-nucleon mg-dependent 
operator with coefficient D2- Until processes involving external pions are computed, the 
contribution of this operator cannot be isolated. Therefore, for demonstrative purposes 
only, we do not modify the depth of the square wells evaluated for physical values of 
in our chiral extrapolation. 

The behavior of the deuteron binding energy as a function of when both the implicit 
and explicit dependences are included is shown in fig. |l^. In this scenario the deuteron 
binding energy decreases as the chiral limit is approached, and the deuteron is unbound 
for ^ 100 MeV. This is precisely the opposite m^-dependence to that found when only 
the explicit dependence is retained. Therefore, at this point in time a reliable calculation 
of the mjr-dependence of the deuteron binding energy cannot be performed. It should 
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Figure 10: The deuteron binding energy as a function of the pion mass that enters both 
explicitly and implicitly in the nuclear potential. For a radius of R = 0.25 fm, an energy- 
independent square-well depth of MVq = 1.03 x 10'^ MeV^ and a square well that depends 
linearly on energy of MV2k^ = 15.25/c^ reproduce the scattering length and effective range 
in the ^Si channel for the physical value of ttIt^ = 139 MeV. 



come as no surprise that the fine-tuning in this channel is significantly disturbed by what 
would appear to be higher-order contributions. Their impact is strongly enhanced by the 
large cancellation between the kinetic and potential energy of the near-threshold bound 
state. 

The ^5*1 — ^Di phase-shifts and mixing parameter are shown in fig. [TTjfor this scenario. 
As required, the ^Si phase-shift is very different in this chiral limit, due to the fact that 




Figure 11: The phase-shifts 6o, 62 and the mixing parameter e\ of the '^S\ — ^Dx coupled 
channels, generated by the parameter set as fig. |T^. The solid curves correspond to the 
Nijmegen phase-shift analysis j{r^, while the two sets of dotted curves correspond to our 
calculated result for the physical value of the pion mass, denoted by in the plots, and 
the chiral limit, as denoted in the plots. 

there is no bound state. On the other hand, the ^Di phase-shift is very similar to that 
computed with only explicit m^-dependence retained. Clearly, the mg-dependence of the 
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OPE couplings makes a perturbatively small contribution to this scattering amplitude, as 
one might expect (there is no pole in this channel). 



4.3 Is an Expansion about the Chiral Limit Sensible? 

We have made the case that an EFT can be developed by performing an expansion about 
the chiral limit. Such an expansion overcomes the formal problem with W power counting, 
that KSW power counting successfully resolved in the ^Sq channel. However, given that 
the deuteron may be unbound in the chiral limit, it is natural to ask whether such an 
expansion makes any sense. The existence or non-existence of the bound state depends 
upon subleading contributions to the scattering amplitude. That is to say, small changes 
in couplings have large effects due to the inherent fine-tuning in the s-wave channels. 

Therefore, in the ^5*1 — coupled channels it may be consistent to perform a chiral 
expansion of observables about a deuteron wavefunction defined at high orders in the 
expansion, as has long been advocated by several groups (see Ref. P, ^) who have 
shown that many observables are relatively insensitive to the choice of regularization. 
For deuteron observables, the EFT would then become an expansion of the interaction 
kernel, up to the order at which the wavefunction is consistently defined. Order by order, 
observables would be independent of i? as i? — if the calculation is performed correctly. 



4.4 Matching to Numerical Results from Lattice-QCD 

In order to match to any foreseeable lattice-QCD results, an extrapolation in rriq will be 
required. What we have presented in this work is the framework required to match to a 
fully unquenched lattice computation, and we have shown that there is currently no control 
of the extrapolation away from the physical value of m^r, due to the fine-tunings that exist 
in nature. However, future lattice computations are likely to be partially-quenched 



as 



opposed to unquenched. Therefore, a partially-quenched EFT will need to be constructed 
in order to correctly perform the mg-extrapolation from the lattice masses, to the physical 
masses. Given our relatively primitive understanding of the EFT for QCD, and also 
considering that the partially-quenched EFT required to perform extrapolations in the 
single-nucleon sector does not exist at this point in time, a significant theoretical effort is 
required to put in place a framework that will enable future lattice-QCD computations 
to make rigorous statements about multi-nucleon systems. 

It is exciting to note that attempts have been made to compute NN scattering param- 
eters with quenched lattice-QCD. Ref. BH] computes the ^5*0 and ^5*1 scattering lengths 



in quenched QCD at /3 = 5.7 with Wilson quarks using Liischer's finite- volume algo- 
rithm, which expresses the energy of a two-particle state as a perturbative expansion 



in the scattering length divided by the size of the box ||39|. Practical computations re- 
quire a sufficiently small box which in turn should be larger than the scattering length. 
Hence Liischer's method is ideal for systems with "natural" scattering lengths, that is 
with size of order a characteristic physical length scale. Of course in the NN system, the 
scattering lengths are much larger than the characteristic physical length scale given by 
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^QCD- Nevertheless one might suppose that away from the physical value of the pion 
mass, the scattering lengths relax to natural values, thus allowing their determination 
from the lattice. In fig. |12| we exhibit the '^5'i scattering length as a function of the pion 
mass as compared to the lattice data from Ref. The left panel has only negligible 



cutoff dependence, by construction. The cutoff dependence exhibited in the right panel 
demonstrates the incomplete nature of the calculation with implicit pion mass dependence 
(the absence of D2). We truncate our curves at ~ 400 MeV as higher-order effects are 
increasingly important in this region. 




200 400 600 200 400 600 

(MeV) (MeV) 

Figure 12: The ^Si scattering length as a function of m.^. In the left panel, the quantities 
qa, M and F.^ are set equal to their physical values for all m,r- I'n the right panel, the 
LO ttLt^ dependence of qa, M and F^^ is included. The dashed (solid) line corresponds 
to a radius of R = 0.25 fm, (0A5fm) and an energy-independent square-well depth of 
MVo = -1.04 X 10"^ MeV^ (-7.82 x lO^MeV^j which reproduce the scattering length for 
the physical value of = 139 MeV. The dots are lattice QCD data, as explained in the 
text. 



5 Conclusions and Thoughts 

We have addressed several important issues concerning the development of a consistent 
and converging EFT to describe multi-nucleon systems. Chiral perturbation theory in 
the single-nucleon sector is a special case of a more general ordering of operators. In the 
presence of more than one nucleon, the momentum expansion is nonperturbative at some 
level to accommodate the fine-tuned scales, while the mg-expansion remains perturbative. 
In other words, the perturbative expansion of nuclear forces is an expansion about the 
chiral limit. In the ^5*0 channel only local operators are treated nonperturbatively, whereas 
in the ^Si — ^Di coupled channels it is necessary to resum the non-local, singular part of 
OPE which survives in the chiral limit. 
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The amount of chiral-limit physics which must be resummed is determined phenomeno- 
logically channel-by-channel. One can always choose to resum the entire momentum 
expansion since perturbative physics remains so even when treated nonperturbatively. 
KSW power counting is the minimal power-counting scheme to describe systems with 
unnaturally-large scattering lengths; i.e. one operator in the momentum expansion is 
summed to all orders. If the effective range is also unnaturally large, it too should be 
resummed and this in turn constitutes a legitimate ordering of operators. W power count- 
ing is a prescription which sums a kernel computed to a given order in chiral perturbation 
theory to all orders. Here there is a subtlety; one must ensure that at a given order in the 
expansion there is a correspondence between divergences and the local operators which 
absorb them. This is obeyed by W power counting in the ^5*1 — ^Di coupled channels but 
not in the ^5*0 channel. A bizarre feature that emerges from the nonperturbative nature 
of the ^5*1 — ^Di coupled channels is the scaling of operators with fractional powers of the 
chiral symmetry-breaking scale, e.g. D2 ~ ^x^^^ 

Our partial higher-order calculations in the '^5'i — '^Di coupled channels suggest that 
an expansion about the chiral limit will converge. However, a full NLO calculation is 
required in order to make a more definite conclusion and to give meaningful predictions 
for the deuteron binding energy in the chiral limit. The use of the square well as a 
short- distance regulator has proved valuable in giving analytic formulas for the RG run- 
ning of the counterterms. However, the necessity of computing processes with external 
gauge fields suggests use of a regulator that manifestly respects gauge invariance, like 
dimensional regularization or Pauli-Villars. An intriguing puzzle remains concerning the 
relationship between square-well regularization and the matching of delta-function in- 
teractions to singular potentials. In taking the limit R ^ one finds that the singular 
potential wavefunctions vanish. This makes it difficult to understand how a delta-function 
interaction can modify the physical asymptotic phase. This puzzle is compounded in the 
coupled-channel system. 

An EFT with dynamical pions is required to describe both elastic and inelastic processes 
involving momenta greater than ~ m^r. For very-low momentum processes -C it is 
appropriate and consistent to work with EFT(7f). All quark-mass dependence is implicit 
in the coefficients of local operators as chiral symmetry is not a symmetry of this EFT. 

One may worry that observables traditionally used to probe nuclear forces, like the 
electric-quadrupole moment, will not be well reproduced by our proposed EFT and will 
be cutoff dependent. However contributions from gauge-invariant, four-nucleon, one- 
photon operators that arise at higher order in the EFT will precisely compensate the 
cutoff dependence and when combined with the one-body contributions from the nuclear 
wave function will reproduce these observables [Q. 

Finally, there is a significant overlap between the approach advocated in this paper, 
and a more mechanical approach to the problem. Haxton and collaborators |^ have 
shown that the nuclear Shell Model can be formulated as an effective theory, and its 
convergence can be dramatically improved by resumming the "short-distance" part of 
the NN interaction and then performing perturbation theory with the remaining "long- 
distance" part. The separation into long- and short- distance is defined numerically by 
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the shape of the NN potential channel-by-channel, rather than by underlying hadronic 
structure. Nonetheless this approach recovers the essential parts of the expansion we are 
advocating in this work. 
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A The Standard Toy Model of Nuclear Physics: Three 
Yukawas 

One can gain significant insight into the perturbative structure of the ^5*0 channel and 
^5*1 — ^-Di coupled channels by studying the behavior of a toy theory of nucleons interacting 
in a potential formed from three Yukawa exchanges [|l^. The potential has the form 



V{r) = Vs{r) + K(r) 

Vs{r) = ap - a„ , = -a^, , (23) 



where the masses 



and couplings 



770 MeV , = 500 MeV , = 140 MeV , (24) 



ap = 14.65 , = 7.00 , = 0.075 , (25) 

are chosen, along with a nucleon mass of M = 940 MeV to approximately reproduce 
both the scattering length, and effective range in the ^5*0 channel, a^^ = —23.7 fm and 
r^^ = 2.57 fm. In making contact with OPE in the real world, the coupling a^^ depends 
upon both the pion mass, m^, and the axial coupling constant, qa, oi-k g\'m\. Therefore, 
perturbing in the pion potential is equivalent to perturbing in either g^-, Tn^^ or both. 
However, it is convenient to think about expanding about the chiral limit, i.e. expanding 
in 5V = V{r] m.,^) — V{r] 0) = K-. One significant difference between this toy theory and 
the real world is that M does not receive contributions from "potential pion" loops, and 
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so M = 940 MeV at each order in perturbation theory. This will not be the case in the 
real world, where M depends analytically and non-analytically on the quark masses. 

Solving the Schrodinger equation in perturbation theory is straightforward. The s-wave 
radial wavefunction, U^^{r), satisfies 



Mdr' 



f/fuii _ {V-E)U^' = . 



full 



(26) 



Writing Uly^ = U^y + + where the superscript denotes the number of insertions 
of it follows that 



Uly - {Vs - E)Uly = 



Mdr^ 



(27) 



with boundary conditions at the origin chosen appropriately for each of the ?7|y. The 
s-wave phase-shifts can be found from the wavefunction at any given order in per- 
turbation theory by simply projecting out the incoming and outgoing asymptotic plane 
waves, e.g. see Ref. [|]. This method automatically provides a unitary phase-shift as 
one is expanding exp {i25) which necessarily involves to all orders. We denote this 
expansion as El. One may perform a further expansion of S^y to any given order in 
that is unitary to a given order in g\, and we denote this expansion as E2. 

At LO, corresponding to no insertions of V^, the "short-distance" phase-shifts resulting 
from El and E2 are identical, as shown in fig. O. One sees sizable deviations between 
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Figure 13: The full phase-shift (solid) and the phase-shift resulting from zero-insertions 
ofVn (dashed) in the three- Yukawa toy model. 



the two phase-shifts below |p| ^ 200 MeV. The phase-shifts at higher orders in perturba- 
tion theory are shown in fig. |14l It is clear that both expansions converge rapidly from 



IpI > 100 MeV. Further, expansion El converges rapidly for all |p|, while E2 converges 
somewhat more slowly in the very low |p| region. 
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Figure 14: The left panel shows the perturbative phase-shifts in the three-Yukawa toy 
model obtained by solving the Schrddinger equation, and using expansion El, while the 
right panel results from expansion E2. The solid curve is the full phase-shift, while the 
dashed, dot-dashed and dotted curves correspond to LO, NLO, and ISPLO, respectively. In 
both cases the ISPLO result essentially lies on top of the full result. 



While we take some joy in demonstrating the perturbative nature of this theory, it 
is still the case that we do not know the underlying short- distance physics responsible 
for the real-world NN interaction and somehow it must be extracted from data. Given 
that we wish to keep the perturbative aspect of the problem, we need to most efficiently 
describe the short-distance part of the interaction, whose exact nature we care little about 
in the low-momentum region. The form of this parameterization has been the source of 
much debate, i.e. W, KSW or some hybrid. However, for this particular problem we 
know exactly what we can do to parameterize the short-distance part only, and that 
is to use effective range theory (ER), which should be valid up to momenta of order 
|p| ~ vti^j'l ~ 250 MeV. Therefore, we fit a polynomial of the form 

|p|cot5, = + + rJIPl' + - , (28) 

a* 2 



to the "short-distance" phase-shift shown in fig. 13, to find 



a' = -3.38 fm , = 2.60 fm , = 0.313 fm , = 0.156 fm . (29) 

A comparison between the "short-distance" phase-shift, and successively better ER ap- 
proximations are shown in fig. |T5|. Retaining only the leading term, i.e. |p| cot 5s = — ^ is 
clearly not such a good approximation for momenta |p| > 40 MeV, and corresponds to the 
underlying perturbative scheme of KSW. However, retaining the short-distance effective 
range contribution is a significant improvement, but still appears to begin deviating from 
the true phase-shift for |p| > 150 MeV. However, this is in no way related to m^r, and is a 
numerical accident. It is this approximation that corresponds to the counting suggested 
in |]T6|, and fleshed-out in Ref. [0. As one proceeds to higher orders in the ER expansion. 



the phase shift is well recovered up to |p| ^ 250 MeV, as expected. 
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Figure 15: Successive effective-range approximations to the short-distance phase-shift. The 
solid curve is the exact short- distance phase-shift while the dashed, dot-dashed, dotted and 
dot-dot- dashed curves correspond to the effective-range expansion including terms up to 
I/qs, Tq, rl and respectively. 




Figure 16: Phase-shifts in the EFT of the three-Yukawa toy model. The dibaryon field 
describing the short- distance physics is expanded out to order |p|^. The left panel shows 
the phase-shift up to one insertion of K-. The solid curve is the full phase-shift at this 
order obtained from the Schrddinger equation, while the dashed, dot-dashed and dotted 
curves correspond to one perturbative insertion of \4 iiT- the EFT expansion including the 
order g\ counterterms up to |pp, |p|^ and |p|^, respectively. The right panel shows the 
phase-shift up to two insertions o/K-- The solid curve is the full phase-shift at this order 
obtained from the Schrddinger equation, while the dashed, dot-dashed and dotted curves 
correspond to two perturbative insertions of in the EFT expansion including the order 
g\ counterterms up to |pp, |p|^ and |p|^, respectively. The g\ counterterms determined 
in the left panel are used in the O {g^ calculation shown in the right panel. 



From the EFT point of view, we can use a generalized dibaryon field to include an 
arbitrary number of terms in the short-distance physics via |p| cot^^- Further, the long- 
distance amplitudes have been computed out to N^LO by Rupak and Shoresh How- 
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ever, there are additional parameters that need to be determined by numerical fitting. 
They occur because the EFT must be matched to the full theory. The short- distance 
behavior of graphs involving all three particles, the tt, a and p, does not allow for a clean 
separation into a dibaryon field defined by the short-distance parameters of eq. ( P^ and 
pion exchange. However, one can perform the separation into a dibaryon field defined by 
the short- distance parameters of eq. (^9|), pion exchange and local finite counterterms of 
the form g'^^ (m,r/'^o-)" -.q'^a (IpI/''^o-)'^ ■> ■■■■ The LO amplitude in the EFT consists solely 
of a dibaryon field with the parameters of eq. (^9|) . The NLO amplitude in the EFT 
consists of the dibaryon field, a single insertion of V-,^ and both momentum-independent 
and momentum-dependent counterterms. As expected, the fit improves as the number 
of counterterms is increased. It is important to note that the size of the counterterms 
is determined by ma- However, if the insertion is not included and the counterterms 
are refit, one finds that their size is now set by m^, and the overall quality of the fit is 
noticeably degraded. At N^LO the pattern continues with the same dibaryon field, two 
insertions of Kr, and more counterterms. The phase-shifts at various orders in VJr are 
shown in fig. |16. 
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